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ABSTRACT 

We present two-dimensional hydrodynamical simulations of both jet- and wind- 
driven models for protostellar outflows in order to make detailed comparisons to the 
kinematics of observed molecular outflows. The simulations are performed with the 
ZEUS-2D hydrodynamical code using a simplified equation of state, simplified cooling 
and no external heating, and no self-gravity. 

In simulations of steady jets, swept-up ambient gas forms a thin shell that can be 
identified as a molecular outflow. We find a simple ballistic bow-shock model is able 
to reproduce the structure and transverse velocity of the shell. Position-velocity (PV) 
diagrams for the shell cut along the outflow axis show a convex spur structure with 
the highest velocity at the bow tip, and low-velocity red and blue components at any 
viewing angle. The power-law index of the mass-velocity relationship ranges from 1.5 to 
3.5, depending strongly on the inclination. If the jet is time-variable, the PV diagrams 
show multiple convex spur structures and the power-law index becomes smaller than 
the steady jet simulation. 

In simulations of isothermal steady wide-angle winds, swept-up ambient gas forms a 
thin shell which at early stages has a similar shape to the shell in the jet-driven model; 
it becomes broader at later times. We find the structure and kinematics of the shell is 
well described by a momentum-conserving model similar to that of Shu et al. (1991). In 
contrast to the results from jet simulations, the PV diagrams for the shell cut along the 
outflow axis show a lobe structure tilted with source inclination, with components that 
are primarily either red or blue unless the inclination is nearly in the plane of sky. The 
power-law index of the mass-velocity relationship ranges from 1.3 to 1.8. If the wind is 
time- variable, the PV diagrams also show multiple structures, and the power-law index 
becomes smaller than the steady wind simulation. 

Comparing the different simulations with observations, we find that some outflows, 
e.g., HH 212, show features consistent with the jet-driven model, while others, e.g., 
VLA 05487, are consistent with the wind-driven model. 

Subject headings: stars: formation — ISM: jets and outflows. 
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1. Introduction 

Protostellar jets and molecular outflows are often associated with the same young stellar objects 
(Lada 1985; Bachiller 1996). They provide unique information about the mass- loss properties of 
young stars and give insight into the star forming process itself. However, the physics connecting 
the two phenomena remains unclear. 

Protostellar jets are generally detected as optical jets in [SII] AA6716, 6731 and Ha (Ray, et 
al. 1996; Reipurth et al. 1999), as infrared jets in H2 (Coppin, Davis, & Micono 1998; Zinnecker, 
McCaughrean, & Rayner 1998), or as radio jets in the radio continuum (Anglada 1995). They are 
also observed in some molecules, e.g., SiO (L 1448, Dutrey, Guilloteau & Bachiller 1997) and CO 
(HH 211, Gueth & Guilloteau 1999). Generally, the jets are highly collimated structures, consisting 
of a series of knots and bow-shocks, extending as far as a few pc out from the central stars (e.g., 
HH 111, Reipurth, Bally, & Devine 1997). 

Molecular outflows are generally seen in CO and some other molecules, as bipolar quasi-conical 
or -parabolic structures surrounding the jet. The observational properties of molecular outflows 
are summarized in Masson & Chernin (1993), Bachiller (1996), Cabrit, Raga & Gueth (1997), and 
Richer et al. (2000). Molecular outflows are much more massive and normally less well collimated 
than the jets. They have a simple power-law mass distribution with velocity (Masson & Chernin 
1992; Stahler 1994). In many cases, the velocities in the molecular outflows increase with distance 
from the source. At low velocities, the outflow appears as a limb-brightened shell surrounding a 
cavity, while at high velocities the outflow sometimes becomes more jet-like. There is also evidence 
that outflow velocities are directed primarily along the major axis of the outflow, with little overlap 
of blueshifted and redshifted emission within the same flow lobe (Meyers- Rice & Lada 1991; Lada &: 
Fich 1996). In a small sample of outflows known to be associated with jets, position- velocity (PV) 
diagrams constructed in a cut along the outflow axis reveal two different kinematic structures: a 
parabolic structure originating at the driving source, and a convex structure with the high velocity 
tip near the H2 bow shock structure (Lee et al. 2000). 

A number of models have been proposed to explain the connection between protostellar jets and 
molecular outflows. Currently, models in which the outflows are driven either by a highly collimated 
jet (Raga & Cabrit 1993; Masson & Chernin 1993; Chernin et al. 1994) or by a wider-angle wind 
(Shu et al. 1991; Li Sz Shu 1996b; Shu et al. 2000) are the most popular. In the former, the bow 
shock originating at the head of the jet sweeps the ambient material into a thin shell which might 
be identified as the molecular outflow. A number of simulations have been carried out to investigate 
the hydrodynamics of protostellar jets (Blondin, Konigl & Fryxell 1989; Blondin, Fryxell & Konigl 
1990; Stone & Norman 1993a,b, 1994; Biro & Raga 1994; Suttner et al. 1997; Smith, Suttner, & 
Yorke 1997; Downes & Ray 1999). Such simulations show that steady, collimated jets are able to 
produce qualitatively the structure and kinematics of some molecular outflows, e.g., NGC 2264G 
(Smith, Suttner, & Yorke 1997). In the wind-driven shell model, ambient material is swept up by 
a radial wide-angle magnetized wind. If the wind has an axial density gradient, the core of the 
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wind can take on the appearance of a collimated jet (Shu et al. 1995; Ostriker 1997; Shang et al. 
1998). If thermal pressure is unimportant, the molecular outflow is momentum-driven (Shu et al. 
1991; Masson k Chernin 1992; Li & Shu 1996b; Matzner & McKee 1999); the resulting structure 
and kinematics are found to be able to describe e.g., HH 111 CO outflow (Nagar et al. 1997). If 
there is little or no cooling of the shocked wind and ambient gas, the dynamics of the outflow can 
be different from the momentum-driven shell model (Frank & Mellema 1996; Delamarter, Frank & 
Hartmann 2000). 

In this paper, we present comparisons of two-dimensional hydrodynamical simulations of both 
jet- and wind-driven outflow models with observed molecular outflows in order to test which aspects 
of each model provides a good fit to the observed kinematics. We also present the first quantitative 
comparison of the structure and kinematics of the swept-up shells found in the simulations (which 
is identified as a molecular outflow) to new and existing analytical models of jet-driven bow shocks 
and wide-angle winds, respectively. Since we concentrate on the dynamics of the large-scale shell, 
our simulations are performed with the ZEUS-2D hydrodynamical code using approximate cooling 
rates, and moderate numerical resolution. We calculate the PV diagrams as well as mass-velocity 
relationships for the shells and compare them with high-resolution CO observations of a number of 
molecular outflows. Calculations are presented for both steady and time-variable (pulsed) jets and 
winds. Our primary goal is to use the simulations to provide a survey of the dynamics of outflows 
for a wide variety of models (winds versus jets, pulsed versus steady), in order to test which models 
best fit the observations. 



2. Equations and Numerical Method 

The two-dimensional hydrodynamic code, ZEUS 2D (developed by Stone & Norman 1992), is 
used to solve the equations of hydrodynamics, 

^+V-(pv) = (1) 

^ + V • (pvv) = -Vp + /?V$ (2) 

de 

— + V • (ev) = -pV ■ v - n 2 A (3) 

where p, v, p, e, and n are the mass density, velocity, thermal pressure, internal energy density, 
and hydrogen nuclei number density, respectively. Some of our simulations use a stratified ambient 
density p; in this case, <I> is the fixed gravitational potential of the ambient material. We adopt an 
optically thin radiative cooling function A for interstellar gas, with the function at high temperature 
from MacDonald &: Bailey (1981) and the function at low temperature from Dalgarno & McCray 
(1972). Since we focus on the dynamics, we assume an equilibrium cooling rate in the simulations, 
in which the ionization fraction of hydrogen is calculated by equating the ionization rate with the 
recombination rate of hydrogen. Simulations with nonequilibrium cooling rate have been performed 
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by e.g., Stone & Norman (1993a,b). Molecular cooling is not included. In our simulations, helium 
is included as a neutral component with n(He) = 0.1 • n, so that n = p/(1.4 • mjj), where uih is 
the mass of atomic hydrogen. An ideal gas equation of state P = (7 — l)e with 7 = 5/3 is used for 
the thermal pressure. Details of the numerical grid and boundary conditions used for the jet and 
wind models are given in §§ 3.1 and 4.1, respectively. 

3. Jet Model 

3.1. Numerical Simulation Parameters 

In our standard model, presented herein, we introduce a perfectly collimated jet into a cold and 
uniform cloud with a density of p = 1.6 x 10~ 20 g cm" 3 (i.e., n = 7 x 10 3 cm -3 ) and a temperature 
of 30 K. The jet is circular with a radius of Rj = 2.5 xlO 15 cm (the typical radius of observed 
optical HH jets; see Ray, et al. 1996), a temperature of 270 K, and a density of pj = 1.6 x 10~ 20 
g cm" 3 (i.e., n = 7 x 10 3 cm" 3 , see Gueth & Guilloteau 1999). Consequently, the ratio of the jet 
density to ambient density, rj, is 1. Speeds in observed molecular jets range from 100 to 500 km 
s _1 , derived from the proper motion of H2 1-0 S(l) knots (Micono et al 1998, Coppin et al. 1998). 
We take the lower limit and adopt Vj = 120 km s _1 for the jet speed, resulting in a mass loss rate 
of 1.2xl0" 7 M yr" 1 . 

The simulations are performed in cylindrical coordinates. We use a computational domain of 
dimensions (R, z) = (1.25 x 10 16 , 1.25 x 10 17 ) cm and a uniform grid of 60 x 600 zones, giving a 
resolution of 2.08 x 10 14 cm, or 12 grid zones per jet radius. Reflecting boundary conditions are 
used along the inner R and inner z boundaries, while outflow boundary conditions are used along 
the outer R and outer z boundaries. For R < Rj along the inner z boundary, inflow boundary 
conditions are used to introduce the jet into the ambient medium. 

We have performed a large number of other simulations in which we vary rj, Rj, and Vj, and 
find results which are qualitatively quite similar to those presented here. 

3.2. Simulation Results of a Steady Jet 

Figure 1 shows the number density, temperature and pressure in the simulation at 650 years, 
with vectors showing the velocity structure. The solid line on the density distribution is calculated 
from a ballistic bow shock model described in a companion paper (Ostriker et al. 2000) and is 
discussed in the next section. With rj = 1, the bow shock speed is about 60 km s _1 . Therefore, the 
jet travels a distance of 1.21 xlO 17 cm in 650 years. Two shocks, a jet shock and a bow shock, are 
clearly seen at the head of the jet in both the temperature and pressure distributions. The bow shock 
interacts with the ambient material and forms a thin shell of shocked ambient material surrounding 
a high-temperature cocoon (consisting of shocked jet gas) which in turn surrounds the jet. The 
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temperature of the shell decreases from 10,000 K at the jet head to about 100 K near the central 
source. In this simulation, only simplified cooling is included. In fact, several other processes, 
such as molecular cooling, collisional ionization of hydrogen, recombination and dissociation of 
molecular hydrogen, and formation of molecular hydrogen on the surfaces of dust grains, would 
also affect the gas temperature determination in the shell. Simulations including these processes 
have been presented by Suttner et al. (1997). In their simulations, there is a significant fraction of 
molecular hydrogen in the shell. Therefore, if molecular cooling were included in our simulation, 
the temperature of the shell would actually be lower. Regions of the shell with temperature lower 
than a few hundred K can be identified with the observed CO molecular outflow, whereas regions 
of the shell with temperature above 1,000 K will produce the observed H2 bow shock. 

Hydrodynamical simulations have been presented by many authors (e.g. Blondin, Konigl & 
Fryxell 1989; Blondin, Fryxell & Konigl 1990; de Gouveia dal Pino & Benz 1993; Stone & Norman 
1993a,b, 1994; Biro & Raga 1994; Suttner et al. 1997; Smith, Suttner, & Yorke 1997; Downes & 
Ray 1999) using different cooling functions and different resolutions. Some of the simulations are 
at higher resolution than this work and thus show much more small scale structure. However, the 
overall shape and size of the shell structure in higher resolution simulations presented elsewhere, and 
also computed by ourselves, is similar to the results we present here, indicating that the underlying 
kinematics is not much different. We will further discuss the effects on the shell due to different 
numerical resolution later in this section. The goal of this paper is to study the overall large- 
scale kinematics of the shell in comparison to CO observations. Detailed high-resolution studies 
of propagating jets necessary for understanding small-scale features in H2 and optical observatoins 
are presented in the papers referred to above. 

The velocity of the material in the shell is almost perpendicular to the shell surface. However, 
the velocity of the material in the cocoon is mostly parallel to the jet axis, indicating that the 
material in the cocoon mostly comes from the jet. The shell structure is corrugated because of 
the variations in the diameter of the jet as it propagates, induced by interaction with shocked 
material in the cocoon (e.g. Blondin, Fryxell k, Konigl 1990). At higher resolution, the shell is 
more corrugated but the overall structure is still the same (see the discussion in next paragraph). 

Figure 2 shows the number density, temperature and pressure along an axial cut through the 
head of the jet. The jet shock and bow shock are clearly seen in the temperature and pressure 
distributions. The shocked ambient material and jet material are compressed by the shocks so that 
the number density peaks between the two shocks. The maximum temperature is 10 4 K at the 
shock fronts. The immediate postshock temperature is T s = 3rhVg /16k, where fh is the average 
mass per particle, or about 10 5 K for 7 = 5/3. For a strong shock, the postshock speed is about 
one- fourth of the shock speed (for 7 = 5/3). The immediate postshock cooling length can be 
approximated as (see Blondin, Fryxell & Konigl 1990), 



where v s and n s are the shock velocity and postshock number density. Since the density increases 
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rapidly toward the interface (contact discontinuity) of the two shocked layers, the cooling rate 
actually increases rapidly toward the interface. For a postshock number density from 10 to 10 6 
cm -3 , the cooling length due to ionic cooling varies between 10 13 to 10 11 cm. This length is 
not resolved in our simulation, therefore the postshock temperature drops to 10 4 K immediately 
behind the shocks. The cooling length behind the shocks in the simulation is thus constrained by 
the numerical resolution, giving a length of about 10 15 cm. Consequently, the shock thickness and 
thus the shell thickness decrease with increasing resolution. In order to investigate the effect of 
numerical resolution of the dynamics of the swept-up shell which surrounds the jet cocoon, we have 
performed a series of simulations with resolution from 4 x 10 14 to 5 x 10 13 cm. We find the shape 
and kinematics of the shell is controlled primarily by the flux of transverse momentum associated 
with dense gas ejected radially from the shocked jet and ambient gas at the head of the jet. We 
find the density of the material ejected from the jet head increases with increasing resolution, so 
that the transverse momentum flux from the jet head does not change significantly with resolution. 
Thus, even though the shell thickness decreases with increasing resolution, the shell structure and 
kinematics do not depend on the resolution. Therefore, our simulation can be used for studying 
the shell structure and kinematics of the jet-driven outflows. 



3.3. Comparison between the Simulation and a Ballistic Bow Shock Model 

In this section, we compare our simulation with the ballistic bow shock model developed 
in a companion paper (Ostriker et al. 2000). In the model, the shell structure and kinematics are 
determined by the transverse momentum flux from the working surface surrounding the hot shocked 
material at the jet head. The transverse momentum flux can be expressed as the ambient mass flux 
flowing into the working surface times a factor proportional to the sound speed c s behind shocks, 
i.e., P r = f37rRjv s pc s , where j3 is a constant and p is the undisturbed ambient density (Ostriker 
et al. 2000). In the simulation, the temperature drops to 10 4 K immediately behind the shocks, 
giving a sound speed of 8 km s _1 . Direct measurement of the transverse momentum flux at R = Rj 
in the simulation gives j3 = 2.2. The value of (3 increases with distance from the working surface 
due to the finite pressure gradient in the shell; also the effective radius of the working surface is 
slightly greater than the jet radius. A larger value of (3 is therefore needed to fit the simulation. 
As can be seen, the shell structure is reasonably described by the model with (5 = 4.1. As a result, 
the shell in our simulation can be expressed as z oc R 3 (Ostriker et al. 2000). This shell structure 
is consistent with that found in the simulations of Smith, Suttner, & Yorke (1997) and Downes & 
Ray (1999). 

Figure 3 presents a comparison of the velocity structure in the simulation and the ballistic 
bow shock model. Figures 3a and 3b show the velocity fields in a frame of reference moving with 
the bow shock for the simulation and model, respectively. In both panels, the velocity is plotted 
as arrows over a gray-scale image of the pressure distribution taken from the simulation. In the 
simulation (Fig. 3a), the hot shocked material at the jet head flows transversely into the shell. As 
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the material moves along the shell, part of the material flows into the cocoon from the inner surface 
of the shell because of the effect of thermal pressure within the shell. The velocity along the inner 
surface of the shell is smaller than the outer surface, producing a velocity gradient across the shell. 
This gradient is a consequence of the fact that the ambient material being swept-up by the shell 
has zero transverse momentum. The bow shock provides an initial transverse impulse; however, 
the swept-up ambient material does not immediately mix with high velocity shocked jet gas ejected 
from the head of the jet, therefore producing the velocity gradient across the shell as seen in Figure 
3a. In Figure 36, we plot the velocity of the newly swept-up material (open arrows) and the mean 
velocity of the material in the shell (closed arrows) from the bow shock model (see Ostriker et al. 
2000). The newly swept-up material lies along the outer surface of the bow shock. The velocity of 
the newly swept-up material has larger magnitude than the mean velocity of the material in the 
shell, indicating that mixing in the shell is not immediate, giving rise to the velocity gradient in 
the shell as seen in Figure 3a. Figures 3c and 3d show the velocity fields for the simulation and 
model, respectively, in the observer's frame plotted as arrows over a gray-scale image of the density 
distribution from the simulation. In the simulation (Fig. 3c), it is clear that there is a gradient in 
the velocity across the shell between the inner and outer shell surfaces, with the outer shell having 
smaller and less forward velocity than the inner shell. In the model (Fig. 3d), the velocity of the 
newly swept-up material (open arrows) is also smaller and less forward than the mean velocity of 
the material in the shell (closed arrows), consistent with the simulation. 

Figure 4 shows transverse velocity, vr, and longitudinal velocity, v z , of the material in the 
shell in the observer's frame for both the simulation (images) and bow shock model (lines). The 
solid lines are calculated with the mean velocity of the shell material and the dashed lines are 
calculated with the velocity of the newly swept-up material from the bow shock model. Both vr 
and v z drop quickly away from the tip of the bow shock, because momentum in the shell is shared 
with the newly swept-up ambient material. The maximum transverse velocity is about 10 km s _1 , 
comparable to the 8 km s -1 sound speed at 10 4 K as predicted by the ballistic bow shock model. 
The transverse velocity of the outer-shell swept-up ambient material is nearly the same as the mean 
transverse velocity of the material already in the shell. Therefore, except near the tip, the dashed 
and solid lines both match the transverse velocity of the simulation reasonably well. As for the 
longitudinal velocity, the solid line (representing mean v z ) matches the high velocity better than the 
dashed line, while the dashed line (representing newly swept-up material) matches the low velocity 
better than the solid line. This is because mixing is more complete near the jet head, so that the 
high velocity is similar to the mean velocity. However, there is almost no mixing in the wing. In 
addition, the material from the working surface partly flows into the cocoon, further reducing the 
mixing. The velocity of the shell material in the wing is thus similar to that expected for newly 
swept-up material. 

Overall, the simple ballistic bow shock model is a good fit for the shell shape and the transverse 
velocity of the shell material. Since the material from the working surface partly flows into the 
cocoon as it moves up the shell and the mixing of the newly swept-up material with material already 
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in the shell is not complete, the model only provides lower and upper limits for the longitudinal 
velocity of the shell material. The lower limit is set by the velocity of the newly swept-up material, 
whereas the upper limit is set by the mean velocity of the shell material. If the numerical resolution 
were much higher so that the shear layer between the swept-up and shell material were resolved, 
mixing might be increased. Moreover, if molecular cooling were included to lower the thermal 
pressure in the shell, there would be less material flowing into the cocoon. In that case, the shell 
material would have a higher forward velocity. Previous simulations which include molecular cooling 
(e.g., Smith, Suttner, k, Yorke 1997; Downes & Ray 1999; Volker et al. 1999) indeed show the 
shell material that has higher forward velocity than that in our simulation. 

3.4. PV Diagrams and MV Relationship 

In our simulation, only a simplified treatment of radiation cooling is included, and moreover, 
the cooling length in the shell is often not resolved. Therefore, the temperature of the shell material 
can not be used to calculate line emission. Instead, our kinematic diagnostics are based on the 
mass rather than the line emission. Since we identify the shell in the simulation as the molecular 
outflow (Lee et al. 2000), we only focus on the shell kinematics and thus exclude the jet and cocoon 
material from our calculations. Since the density contrast between the shell and cocoon material 
is large (see Figure 1), we can define a boundary between the shell and cocoon material and mask 
out the cocoon and jet material. Unavoidably, a little cocoon material close to the shell will be 
included in our calculations; however, the mass contribution of that cocoon material is small. 

Figure 5 shows position-velocity (PV) diagrams for the shell in the simulation cut along the 
outflow axis at three inclinations (i) to the plane of the sky. The jet and cocoon material are 
masked and excluded from the calculations. The PV diagrams are calculated by projecting the 2- 
dimensional density distribution into 3-dimensions and summing the mass along each line of sight 
for any given velocity. The solid lines are PV diagrams calculated from the ballistic bow shock 
model using the mean velocity of the shell material, while the dashed lines using the velocity of 
the newly swept-up ambient material. It is clear that the shell material in the wings is dominated 
by the unmixed shocked ambient material, so that the observed velocity is similar to the observed 
velocity of the newly swept-up ambient material. The PV diagrams are similar to that seen in other 
simulations: they consist of convex spur structures with the highest velocity at the tip (Smith, 
Suttner, & Yorke 1997; Downes & Ray 1999). At i = 0°, where the outflow axis is in the plane of 
the sky, the diagram is symmetric about zero velocity. As the inclination increases (with the jet 
pointing away from the observer), the blueshifted shell material and low- velocity redshifted shell 
material moves toward zero velocity, while the high-velocity redshifted shell material moves away 
from zero velocity. Therefore at i = 60°, the blueshifted shell material and low- velocity component 
of redshifted shell material merge with the ambient material, leaving the high- velocity shell material 
mostly visible on the redshifted side. 

Figure 6 shows mass-velocity (MV) relationship for the shell material at the same three in- 
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clinations. The MV relationship is computed by summing the mass for any given velocity. Both 
the redshifted (open square) and blueshifted (filled square) masses are shown. The curves in the 
figure are not smooth since the shell itself is not smooth. The dashed lines are the fits to the 
redshifted mass with a power-law mass velocity relationship, dM/dv b s = m(v t, s ) oc , where 
the power-law index, 7, is indicated at the upper right corner in each panel. As found by Downes 
k. Ray (1999), our result shows that 7 depends strongly on the inclination; 7 is 3.44 at % = 0° 
and decreases to 2.01 at i = 60°. 7 is also consistent with that found by Smith, Suttner, & Yorke 
(1997) if their channel line intensity scales with dm(v) rather than v 2 dm(v). Notice that what is 
actually observed is a variation in CO line intensity with velocity. If the CO line is optically thin 
and the temperature of the gas is constant and higher than the excitation temperature of the line, 
CO line intensity is directly proportional to mass (see, e.g. McKee et al. 1982). Therefore, our 
mass-velocity relationship can only be applied to the observations if the CO line is optically thin 
and the temperature of the shell is constant. Readers should consult e.g. Smith, Suttner, & Yorke 
(1997) for a more sophisticated calculation of CO line intensity-velocity relationship. 

Since the bow shock model (Ostriker et al. 2000) can reasonably reproduce the transverse 
velocity in the simulation, the model can also account for the MV relationship at i = 0, where the 
observed velocity depends on the transverse velocity only. The solid line at i = 0° in Figure 6 is 
calculated from the bow shock model and matches the MV relationship of the simulation well. 



3.5. Comparison with a Pulsed Jet Simulation 

Optical and infrared observations reveal that protostellar jets usually consist of a series of well 
aligned knots and internal bow shocks moving away from the source at high velocity (Coppin, Davis, 
& Micono 1998; Micono et al. 1998). These knots and bow shocks are thought to be associated with 
internal working surfaces produced by time variability in the jet (Reipurth & Graham 1988; Raga 
et al. 1990). In this section, we present a simulation of a time- variable (pulsed) jet and compare it 
with the steady jet simulation. The jet velocity and density are assumed to be 

v'j = Vj (l + Asin^) (5) 

p^ = Pj ./(l+^sin^) (6) 

where A is the amplitude of the variation and P is the period. The momentum of the jet is 
constant over the variation. We use A = 0.5 and P = 310 years, allowing a strong variation of the 
jet momentum to demonstrate the effects of the internal working surface. 

Figure 7 shows a gray-scale image of the density distribution in the simulation at 610 years, 
along with PV diagrams and the MV relation at two inclinations. At a time of 310 years, two internal 
working surfaces are formed in the jet, producing a knot and an internal bow shock structure. Due 
to thermal pressure gradients, shocked jet material is ejected radially from each internal working 
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surface into the cocoon. The ejected material initially appears as a knot. As the working surface 
travels down the jet axis, the knot grows into an internal bow shock surface and eventually interacts 
with the shell. Since the inertia of the shell is large, the shell structure is not significantly affected by 
the internal bow shock (see e.g., Stone & Norman 1993b; Biro h Raga 1994; Suttner et al. 1997). 
The internal bow shock, which consists of jet material, has a higher forward velocity than the 
leading bow shock surface, because it interacts only with the low-density, forward-moving cocoon 
and not the dense, stationary ambient gas. From the PV diagram at i = 0°, it is clear that the 
velocity of the internal bow shock surface decreases slowly with R before it interacts with the shell, 
again because there is not much material in the cocoon. As the bow shock surface interacts with 
the shell, momentum is transfered to the shell material and the velocity of the shell increases. As 
a result, the velocity structure of the shell is broken at the interaction point into two convex spur 
PV structures, one associated with the leading bow shock and one with the internal bow shock. 
The internal bow shock surface has a higher forward velocity and appears as high velocity gas in 
the PV diagram at i = 30°. The mass at high velocity in this simulation is higher than that in the 
steady jet simulation. Therefore, the mass-velocity relationship in the pulsed jet simulation has a 
smaller power-law index than that in the steady jet simulation (compare Figures 6 and 7). 



3.6. Steady Jet in a Highly Stratified Ambient Medium 

We have also computed a steady jet in a highly stratified ambient medium with density dis- 
tribution (in cylindrical coordinates) 

where p c is the density at z = and z c is the effective core-radius of the ambient medium. We use 
p c = 1.6 x 1CT 19 g cm~ 3 (i.e., n = 7 x 10 4 cm~ 3 ) and z c =1.25 xlO 17 cm. The jet is underdense 
(77 = 0.1) at z = with r\ increasing as the jet propagates down the jet axis. Due to the strong 
cooling of the shocked material at the head of the jet and the cocoon material, we find the structure 
and kinematics of the swept-up shell in an underdense jet are more or less ballistic and similar to 
that of the steady jet in an uniform ambient material with 77 = 1 (see §3.2). Thus, stratification of 
the ambient medium above does not modify the kinematic signature of the jet-driven outflow. 

4. Wide-angle Wind Model 

4.1. Numerical Simulation Parameters 

The wide-angle wind model is an alternative to the bow shock of a collimated jet. In this 
model, a molecular outflow is ambient material swept-up by a radial wind from a young star. The 
wind could be stratified in density if it emerges from an extended disk (e.g. Ostriker 1997) or by 
the action of latitudinal magnetic stresses in a rotating protostellar x-wind (Shu et al. 1995). In 
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a spherical coordinate system (r,0,cj)), Shu et al. (1995) have shown that for an x-wind model the 
wind density can be approximated as 

p w oc l/(rsin6>) 2 . (8) 

This is also true for more general force- free winds (Ostriker 1998; Matzner & McKee 1999). We 
have performed a number of time-dependent hydro dynamical simulations of the evolution of wide- 
angle wind in stratified ambient medium. The wind is introduced as a boundary condition along a 
spherical shell of radius r w , with density distribution 

(9) 



,w r2(sin 2 # + e) 

where p wo is a constant and e is a small value avoiding a singularity of the wind density at the pole 
(0 = 0). In the x-wind model, the wind velocity v w is approximately the same at all angles (Najita 
& Shu 1994). Here, we assumed the wind velocity decreases toward the equator (6 = ir/2) and is 
given by 

v w = v wo cos9 (10) 

where v wo is the velocity at the pole. The mass loss rate of the wind is set to that in the jet 
simulations, thus 

rn/2 

4vr / p w v w rl sm9d6 = 2p j nR 2 j v j = 1.2 x 1O' 7 M yr' 1 (11) 
Jo 

Letting v wo = Vj, r w = Rj and e = 0.01, we have p wo = 3.5x 10' 21 r^, g cm -1 , giving p w = 3.5x 10' 19 
g cm' 3 (i.e., n = 1.5 x 10 5 cm' 3 ) at the pole, and p w = 3.5 x 10' 21 g cm' 3 (i.e., n = 1.5 x 10 3 
cm' 3 ) at the equator. 

For the ambient medium, a flattened torus with density 

P° = ~2 ( 12 ) 

is used, where p ao is a constant. This density distribution is appropriate for magnetized cores 
(see e.g., Li & Shu 1996b). In the simulations, p ao is set to 1.6 x lO' 18 ^ g cm' 1 , which leads 
p a = 1.6 x 10' 18 g cm' 3 (i.e., n = 6.8 x 10 5 cm' 3 ) and rj = 2.2 x 10~ 3 at the equator when r = r w . 
We also added a fixed gravitational potential $ to keep the ambient material stable. This is the 
most favorable distribution of the ambient material for the wind model. Analytically, it has been 
shown to produce not only the correct mass- velocity relationship, but also the collimation required 
for molecular outflows (Li & Shu 1996b; Matzner & McKee 1999). 

In the x-wind model, the wind and ambient material are actually magnetized (Shu et al. 1994, 
1995; Li k, Shu 1996b). For a magnetized gas with a low ionization fraction, a non-dissociative 
C-shock can occur for shock velocities up to 40 — 50 km s" 1 (Hollenbach 1997). Therefore, a C- 
shock interaction is likely to take place in the wings of the outflow lobe, with a J-shock interaction 
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near the tip of the lobe. At temperatures below 10 4 K, the cooling of the shocked material in the 
wings is thus dominated by molecular cooling (Hollenbach 1997). Since molecular cooling is very 
efficient for the density used in our simulation, an isothermal state function is used for the thermal 
pressure. The effect of magnetic fields on the dynamics of wide-angle winds will be presented in a 
future publication. 

The simulations of this model are performed in a spherical coordinate system. We use a 
computational domain of dimensions (r, 0) = (2.5 x 10 15 — 1.25 x 10 17 cm, tt/2) and a uniform 
grid of 600 x 400 zones, giving a resolution of 2.04 x 10 14 cm in r and 4xl0~ 3 radians in 0. 
Reflecting boundary conditions are used along the inner (where = 0) and outer (where 
9 = it/2) boundaries, while outflow boundary conditions are used along the outer r boundary. 
Inflow boundary conditions are used for inner r boundary to introduce the wind into the ambient 
medium. 



Figure 8 presents the density distribution of an isothermal wind on a logarithmic scale at 390 
years, with vectors showing the velocity structure in the shell. The solid line indicates the predicted 
shape from a simple momentum-drive shell model derived in the next section. In the simulation, 
the temperature is set to 100 K, a typical temperature of the CO outflow in the observations (Fukui 
et al. 1993). The central high density "jet" spreads with a small opening angle because we have 
used a small factor, e, to avoid the singularity of the wind density at the pole. The shell is elongated 
in the polar direction, due to the stratification of the ambient material and the variation in angle 
of the wind thrust. The shell everywhere expands in the radial direction. At a temperature of 
100 K, the thermal pressure in the shell is negligible compared to the ram pressure of the wind. 
Consequently, the shocked wind and ambient materials merge into a thin shell of shocked material. 

4.3. Comparison Between the Simulation and A Momentum-driven Shell Model 

Assuming immediate cooling behind shocks and no relative motion of the shocked material, 
Shu et al. (1991) and Li & Shu (1996b) derived the shape and kinematics of the swept-up shell 
by considering mass and momentum conservation in each angular sector. In their calculations, the 
wind mass added to the shell is not included. Following the same approach, we rederive the shape 
and kinematics of the shell including the wind mass added to the shell below. The mass (ambient 
mass + wind mass) and momentum fluxes added to the shell per second per steradian are 



4.2. Simulation Results for an Isothermal Steady Wind 



dM. 



s 




r 2 w p w {y w - v s ) 



(13) 



dt 



d(M s v s ) 
dt 



rlp w (v w - v s )v. 



w 



(14) 
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where M s is the mass flow into the shell and v s is the shell velocity. Since the ambient material 
has a density varying as r -2 , the dilution of the wind ram pressure is exactly compensated as the 
outflow expands away from the source, resulting in a swept-up shell which proceeds outward with 
a constant velocity along any radial line (Shu et al. 1991; Masson &; Chernin 1992). Therefore v s 
is constant and 

dM. 

= r w p w (v w - v s )v w (15) 



Vs dt 

Substituting equation (13) into equation (15) and solving for v s , we have 



1 +Vw ' 



2 

where tj w = The shell structure is then given by 



r s = v s t + r w = — Vwt U9 + r w (17) 

i + v™ 

The maximum width of the shell can also be found at angle sin# = [2 + (■^-) 1 / 2 ] -1 / 2 as 

L v pwo ' J 

W w « W(^) 1/4 (18) 

Pao 

for W w 3> r w and e = . Since the length of the shell is ~ v wo t, the width-to-length ratio is about 
r£wo\i/4: ^ constant with time. The ratio of the total longitudinal momentum flux to the total 

V Pao ' ° 

transverse momentum flux in the shell can be expressed as 

-1 + i±^arctan^ 

n « fi= (19) 

where a = « Ems.. In our simulation, ((has.) 1 / 4 0.22 and 7Z 4.5, therefore, the wind-driven shell 

Y Pao ' V Pao ' 

is elongated in the polar direction. 

The solid line in Figure 8 indicates the shell shape calculated from the momentum-driven shell 
model. As can be seen, the model provides a good fit. Due to finite thermal pressure gradients, the 
shell in the simulation is a little larger than the model. Figure 9 shows the transverse velocity, vr, 
and the longitudinal velocity, v z , along the shell, with solid lines indicating the predicted velocities 
from the model. Again, the model provides a good description of the velocity structure. The 
transverse velocity increases from zero at the source and then decreases to zero at the tip, while 
the longitudinal velocity increases linearly with distance. In contrast to the jet simulation, the 
transverse velocity vanishes at the tip because the thermal pressure is negligible in the isothermal 
simulation at such a low temperature. 

Comparing Figure 8 with Figure 1, we see that the wind-driven shell is wider than the jet- 
driven shell, although the wind and jet both have the same mass loss rate in the simulations. The 
width-to-length ratio for the jet-driven shell is 2(3Pc s /v s ) 1 / 3 \Rj/z\ 2 / 3 and decreases with increasing 
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z (Ostriker et al. 2000). The width-to- length ratio for the wind-driven shell, on the other hand, is 
^£i£o -ji/4 an( ^ i n d e p enc i en t of z. Therefore, even with optimized ambient medium stratification most 
favorable for narrow outflows in the wind model and wide outflows in the jet model, the wind-driven 
shell always becomes wider than the jet-driven shell beyond a certain distance of z; in our case, 
this occurs at about 30 Rj. 



4.4. PV Diagrams and MV Relationship 

Figure 10 shows PV diagrams for the shell material cut along the outflow axis at three inclina- 
tions, with lines indicating the model calculations. The PV diagrams are significantly different from 
that seen in the jet simulations, showing a lobe structure that is tilted with inclination such that 
primarily red or blue, but not both, emission is seen if the inclination is nonzero. At an inclination 
i, the projected distance of the shell material is 

d = r s (cos 9 cos % — sin 9 sin i cos <p) (20) 

and the observed velocity is 

v bs = v s (sm 9 cos (f> cos i + cos 6 sin i) (21) 

where <fr is the azimuthal angle in the shell. Therefore the PV diagram along the major axis is 
given by 

d = ( , v ° bs ^ 1- r w )(cos9cosi — sin 9 sin i cos 4>) (22) 

sm 9 cos <p cos i + cos dsini 

where <p = for the front wall and <p = tt for the back wall of the shell. As can be seen in Figure 
10, the PV diagrams from the simulation are well fit by the model. 

Figure 11 shows the mass- velocity (MV) relationship for the shell material. Both the redshifted 
(open square) and blueshifted (filled square) masses are shown. The dashed lines are the fits to the 
redshifted mass with a power-law mass velocity relationship, m(v b s ) ex v ob J s , where the power-law 
index, 7, is indicated at the upper right corner in each panel. The solid lines are from the model 
calculations. The swept-up mass (ambient mass + wind mass) per steradian, m((i, 0), is 

m(fi,(f)) = j^^-dt = rlp w v w r/- 1/2 t (23) 



where ji = cos 9. Therefore, the mass per unit velocity is 

# 
d/i 

dv bs 

—dfi 



m(v obs ) = J m(fi,(f>) 



= J m{n,4>)- 

(24) 



v s sin 9 sin <p cos i 
-rlp w v w T]w 1/2 tdp 



\J v% sin 2 9 cos 2 i — (v b s — v s cos 9 sin i) 2 
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In the simulation, the wind sweeps up a slightly larger shell than predicted by the model (see Figure 
8), therefore, the mass from the model is multiplied by a factor of 1.3 to match the simulation. 
Again, the model matches the MV relationships very well. The slope of the curve is steeper at 
higher velocity. Fitting to the mass beyond a few km s _1 gives a 7 from 1.3 to 1.8, a little smaller 
than the original x-wind model, which is about 2 (see, e.g., Li & Shu 1996b; Matzner & McKee 
1999). 



4.5. Comparison with an Isothermal Pulsed Wind Simulation 

We have also computed the evolution of an isothermal time-variable (pulsed) wind for com- 
parison to the simulation of a isothermal steady wind. As in the pulsed jet simulation, the wind 
velocity and density are assumed to be 

v' w = v w (l + Asm—) (25) 

27rf 

p' w = Pw /(l + Asm— ) (26) 

where A is the amplitude of the variation and P is the period. Again, the momentum of the wind 
is constant over the variation. In the simulation, A = 0.5 and P = 115 years. 

Figure 12 shows a gray-scale image of the density distribution in the simulation at 296 years, 
along with PV diagrams and the MV relation at two inclinations. As the wind velocity changes, the 
fast wind overtakes the slow wind, producing internal shocks inside the shell. As was seen in the 
pulsed jet simulation, the shell also shields the internal surfaces from directly interacting with the 
ambient material. The polar region of the internal shocks moves down the axis without interacting 
with the shell until it reaches the tip of the shell. The structure of the internal shocks is thus 
determined by the intrinsic angular distribution of the wind velocity, as indicated by the solid lines, 
and is flatter further out from the source. This structure of the internal shocks is much flatter than 
the observed internal H2 bow shocks in e.g., HH 212 (Zinnecker, McCaughrean, & Rayner 1998) 
and HH 111 (Coppin, Davis, & Micono 1998). To account for the observed shape of internal H2 
bow shocks in these outflows, the temporal variation in the wind model would need to be confined 
to a small angle along the outflow axis, and the wind speed would also need to decrease more 
rapidly toward the equator. In the simulations, the wings of the internal surfaces interact with the 
shell and transfer momentum to the shell, producing multiple structures in the PV diagrams, each 
associated with an internal surface. In the PV diagrams, the polar region of the internal surfaces 
appear almost as straight lines. The power-law index of the MV relation is smaller than the steady 
wind simulation, indicating there is more mass at higher velocity. 



5. Comparison with Observations 



Currently, the best example for comparison with the features evident in the jet-driven model 
is probably the HH 212 outflow (Zinnecker, McCaughrean, & Rayner 1998). In Figure 13, the 
left panel shows CO emission contours (Lee et al. 2000) plotted over an H2 image (Zinnecker, 
McCaughrean, & Rayner 1998), and the right panel shows the PV diagram of the CO emission 
cut along the jet axis (Lee et al. 2000). The jet itself is seen as a series of knots and bow shock 
structures in H2 emission along the outflow axis. The CO emission surrounds the H2 emission. 
This morphological relationship can be produced in a jet simulation, in which the H2 emission is 
produced near the bow shock surface and CO emission is produced in the wings and around the 
bow shock surface. The PV diagram of the CO emission in HH 212 along the jet axis shows a 
series of convex spur structures on both redshifted and blueshifted sides with the highest velocity 
near the H2 bow tips, qualitatively similar to the PV diagram at zero inclination in our pulsed jet 
simulation (see Fig. 7). Notice that the PV diagram of the outflow is a little asymmetric due to 
the inclination effect. At i = 0°, the PV diagram in the ballistic bow shock model can be expressed 
as v bs = [/3c s v 2 i? 2 /9] 1//3 |z|~ 2//3 , where \z\ is the projected distance from the bow shock. For an 
internal surface, this relationship becomes v Q b s = [f3c s Av(v s i — v e )R 2 / (\Sp e / pj)] l ^\z\^ 2 ^ , where 
v s i is the internal shock speed, At; is velocity jump across the internal shock, v e and p e are the 
velocity and density of the cocoon material (Ostriker et al. 2000). The solid lines in the Figure 
13 are calculated with [f3c s /\v(v s i — v e )R 2 /(18p e / pj)] = 3650 arcsec 2 (km/s) 3 and Rj in unit of 
arcsec. With Rj = 0.2" (i.e., 100 AU at 460 pc) and [3c s = 32 km/s, we have A "^/ ~" e) = 5.1 x 10 4 
(km/s) 2 . For HH 212, v si is about 75 km s" 1 (Davis et al. 2000). If v si - v e « 50 km s -1 , Av » 50 
km s _1 , then p e /pj = 0.05. The close agreement between the observed properties of the HH 212 jet, 
and the parameter values needed to fit the observations with the jet-driven outflow model, leads 
us to conclude this molecular outflow may be driven by the underlying jet. 

Two other good examples for comparison with the jet-driven model are the HH 240/241 outflow 
(Lee et al. 2000) and the NGC 2264G outflow (Lada & Fich 1996). The CO emission of these 
outflows show bow shock structures. The PV diagrams of the HH 240/241 (the western lobe) and 
NGC 2264G CO outflows along the major axis also show the similar PV structures to that seen 
in our simulations at an inclination of 60°. The actual inclinations of these outflows depend on 
the mixing of the shell material. Smaller inclination is required to reproduce the PV structures of 
these outflows if the mixing is more complete. 

A good example of the wind-driven model is the VLA 05487 outflow (Lee et al. 2000). In 
Figure 14, the left panel shows the CO emission contours (Lee et al. 2000) plotted over the H2 
image (Garnavich et al. 1997), and the right panel shows the PV diagram of the CO emission cut 
along the jet axis (Lee et al. 2000). Since there are other two outflows to the south contaminating 
this outflow, we concentrate on the northern lobe. The northern lobe of the CO outflow forms 
a conical structure around jet-like H2 emission. The PV diagram along the major axis shows a 
parabolic PV structure extending out from the source. The CO emission structure and PV structure 
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have been reasonably modeled with a radially expanding parabolic shell with a Hubble-law velocity 
structure, which is derived from the x-wind model (Lee et al. 2000). This outflow is large, extending 
at least 8' to the north (Reipurth & Olberg 1991), the observations in Lee et al. (2000) only revealed 
the inner part of the CO molecular outflow. Therefore the emission structure and PV structure 
of this outflow both appear parabolic. The solid line in the PV diagram is calculated from our 

l 27 fkni/sl^ 

modified x-wind model with i = 18°, t = 8100 years and ^ = - — ^ — at a distance of 460 pc. 
This gives v = 0.27 km s _1 arcsec -1 and C = 0.24 arcsec -1 using Equation 1 in Lee et al. (2000), 
consistent with their work. As can be seen, the model can account for the PV diagram reasonably 
well. 

In the observations, the power-law index of the mass-velocity relationship, 7, is calculated 
directly from the CO emission of the molecular outflows. Previously, 7 is known to be about 1.8 if 
the emission is optically thin or if constant optical depth emission was assumed for the CO emission 
(Masson & Chernin 1993; Cabrit, Raga & Gueth 1997). This 7 is consistent with the wind-driven 
model. However, if the 12 CO optical depth is corrected using a velocity-dependent fit to the ratio 
r ( 12 COJ = 1 - 0)/r( 13 COJ = 1-0), then 7 is found to be between 2 and 4 (Yu 1999), which is 
more consistent with the jet-driven model. Further observations are needed to clarify the actual 
mass- velocity relationship. 

6. Conclusions 

We have performed a systematic study of the basic properties of both jet-driven and wind- 
driven molecular outflows using time-dependent hydrodynamic simulations and simple analytic 
models. The main conclusions are the following: 

1. The structure and transverse velocity of the swept-up shell in a steady jet propagating into 
a uniform ambient medium can be reasonably reproduced by a ballistic bow shock model, 
presented in a companion paper (Ostriker et al. 2000). Since some material flows into the 
cocoon from the inner shell surface, and since mixing of newly swept-up material with material 
already in the shell is not complete, the ballistic bow shock model only provides lower and 
upper limits for the longitudinal velocity of the shell in the simulation. 

2. The PV diagrams along the outflow axis for the shell material swept-up by a steady jet show 
a convex spur structure with the highest velocity at the bow tip. Low- velocity shell material 
is relatively symmetric in red and bluesides at any inclination; high velocity material is one- 
sided (either red or blue). The power-law index of the mass- velocity relationship ranges from 
1.5 to 3.5, depending strongly on the inclination. 

3. In a pulsed jet, the internal bow shocks can affect the shell kinematics, producing multiple 
convex structures in the PV diagram. The internal bow shocks consist of jet material, and 
have a higher forward velocity than the leading bow shock, producing a high velocity gas 
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signature in the PV diagrams even at small inclination. The leading bow shock structure is 
not strongly affected by internal bow shocks that collide with it obliquely. 

4. The structure and kinematics of the shell in an isothermal steady wide-angle wind simulation 
can be well described by a momentum-driven shell model similar to that of Shu et al. (1991). 
Since the thermal pressure is small compared to the ram pressure of the wind, the shocked 
wind and ambient materials merge into a thin shell of shocked material, so that mass and 
momentum are conserved in each angular sector. 

5. In the steady wind simulation, the PV diagrams cut along the outflow axis show a lobe 
structure tilted with inclination, so that any given lobe is primarily red or blue except if the 
axis is nearly in the plane of the sky. Our wind models can not produce the spur-like features 
seen in the jet simulations and a number of observed systems. The power-law index of the 
mass-velocity relationship ranges from 1.3 to 1.8. 

6. In a pulsed wind, the polar region of the internal surfaces is flatter further out from the source 
and the wings of the internal surfaces interact with the shell producing multiple straight 
structures in the PV diagrams. The internal surface consists of the wind material that has a 
higher forward velocity than the leading surface, producing high velocity gas signature in the 
PV diagrams even at small inclination. 

7. The overall width of jet-driven shell is smaller than that of wind-driven shell, even when the 
ambient medium stratification is most favorable for narrow outflows in wind model and wide 
outflows in jet model. 

There are clear differences between the PV diagrams and the power-law index of the mass- 
velocity relationship between these two models. Comparing to observations, we find that some 
outflows, e.g., HH 212, are consistent with the jet-driven model, while others, e.g., VLA 05487, 
are consistent with the wind-driven model. Since the internal surface in the wind model is much 
flatter than the observed H2 bow shocks, the temporal variation in the wind model would need to 
occur within a small angle along the outflow axis to produce such discrete curved structures, and 
a significant decrease in the wind velocity toward the equator would also be required. 

While none of our simple dynamical models of protostellar outflows are consistent with all of 
the detailed kinematics of individual sources, this work shows that broad similarities with charac- 
teristic signatures are easily identified. The actual jet, wind, and ambient material are undoubtedly 
more complicated than our simple models, due to the intrinsic properties of the systems (e.g., mag- 
netic fields, rotation, density and velocity fluctuations in ambient gas). Such complexities could 
potentially explain why some sources (e.g. HH 111, Lee et al. 2000) evidence signatures of both 
wind and jet driving in their outflow kinematics. It is possible that the differences among observed 
outflows represent an evolutionary effect. Further simulations are needed to investigate these and 
other intriguing questions. 
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Fig. 1. — Distributions of number density, temperature and pressure in logarithmic scale at 650 
years in the steady jet simulation. The vectors show the velocity structure. The solid line plotted 
over the number density is the shell shape calculated from a ballistic bow shock model derived 
by Ostriker et al. (2000). The gray-scale wedges along the right side indicate the values for the 
gray-scale images. The z-axis is along the jet flow axis and there is an assumed circular symmetry 
about the z-axis. 

Fig. 2. — Number density, temperature, and pressure along an axial cut through the head of the 
jet in the steady jet simulation. 

Fig. 3. — Comparisons of the velocity structures of the shell material between the steady jet 
simulation and ballistic bow shock model. Panels a and b show the velocity structures for the 
simulation and model, respectively, in the bow shock frame plotted over the pressure distribution 
of the simulation. Panels c and d show the velocity structures for the simulation and model, 
respectively, in the observer's frame plotted over the density distribution of the simulation. In 
panels b and d, the open arrows are calculated with the velocities of the newly swept-up material 
and closed arrows are calculated with the mean velocities of the shell material from the model. 

Fig. 4. — Comparisons of the transverse velocity and longitudinal velocity of the shell material in the 
steady jet simulation (image) with the ballistic bow shock model (lines). Solid lines are calculated 
with the mean velocity of the shell material. Dashed lines are calculated with the velocity of the 
newly swept-up material. 

Fig. 5. — PV diagrams for the shell material at three inclinations cut along the outflow axis for 
the steady jet simulation, i is the inclination of the outflow to the plane of the sky. Solid lines 
are calculated using the mean velocity of the shell material. Dashed lines are calculated using the 
velocity of the newly swept-up material. Dotted lines indicate the zero velocity. 

Fig. 6. — Mass-velocity relationships at three inclinations for the steady jet simulation. Both the 
redshifted (open square) and blueshifted (filled square) masses are shown. The dashed lines are the 
fits to the redshifted mass with a power-law mass velocity relationship, where the power-law index, 
7, is indicated at the upper right corner in each panel. The solid line at i = 0° is calculated from 
the ballistic bow shock model. 

Fig. 7. — Distribution of number density in logarithmic scale at 610 years in the pulsed jet simu- 
lation, with vectors showing the velocity structure. PV diagrams and MV relationships at i = 0° 
and i = 30° are also shown. The gray-scale wedge along the right side indicates the values for 
the number density. Dashed lines in the MV relationships are fits to the redshifted mass with a 
power-law mass velocity relationship, where the power-law index, 7, is indicated at the upper right 
corner. 

Fig. 8. — Distribution of number density in logarithmic scale at 390 years in the isothermal steady 
wind simulation, with vectors showing the velocity structure in the shell. The solid line is the shell 
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shape calculated from the momentum-driven shell model. The gray-scale wedge along the right 
side indicates the values for the number density. 

Fig. 9. — Transverse velocity, vr, and longitudinal velocity, v z , of the shell material. The solid 
lines indicate the velocity calculated from the momentum-driven shell model. 

Fig. 10. — PV diagrams for the shell material at three inclinations cut along the outflow axis, i is 
the inclination of the outflow to the plane of the sky. The solid lines are the PV diagrams calculated 
from the momentum-driven shell model. 

Fig. 11. — Mass- velocity relationships at three inclinations. Both the redshifted (open square) 
and blueshifted (filled square) masses are shown in the figure. The dashed lines are the fits to 
the redshifted mass with a power-law mass velocity relationship, where the power-law index, 7, 
is indicated at the upper right corner in each panel. The solid lines are calculated from the 
momentum-driven shell model. 

Fig. 12. — Distributions of number density in logarithmic scale at 296 years in the isothermal pulsed 
wind simulation. The vectors show the velocity structure in the shell and internal surfaces. The 
gray-scale wedges along the right side indicate the values for the gray-scale images. PV diagrams 
and MV relationships at i = 0° and i = 30° are also shown. Dashed lines in the MV relationships 
are fits to the redshifted mass with a power-law mass velocity relationship, where the power-law 
index, 7, is indicated at the upper right corner. 

Fig. 13. — CO emission, H2 emission and PV diagram of the HH 212 outflow. On the left is 
the CO emission contours (Lee et al. 2000) plotted over the gray-scale image of the H2 emission 
(Zinnecker, McCaughrean, h Rayner 1998). The outflow has been rotated clockwise by 23°, so that 
the outflow axis is N-S oriented. On the right is the PV diagram of the CO emission cut along the 
jet axis provided by Lee et al. (2000). Vertical dashed line indicates the ambient velocity around 
the outflow. Horizontal dashed line indicates the position of the driving source of the outflow. The 
solid lines plotted over the PV diagram are calculated with the ballistic bow shock model. 

Fig. 14. — CO emission, H2 emission and PV diagram of the VLA 05487 outflow. On the left is 
the CO emission contours (Lee et al. 2000) plotted over the gray-scale image of the H2 emission 
(Garnavich et al. 1997). The outflow has been rotated counterclockwise by 4°, so that the outflow 
axis is N-S oriented. On the right is the PV diagram of the CO emission cut along the jet axis 
provided by Lee et al. (2000). Vertical dashed line indicates the ambient velocity around the 
outflow. Horizontal dashed line indicates the position of the driving source of the outflow. The 
solid line plotted over the PV diagram is calculated with the modified X-wind model. 
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